Calculated NMR T2 relaxation due to vortex vibrations in cuprate superconductors 
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We calculate the rate of transverse relaxation arising from vortex motion in the mixed state of 
YBa2Cu307 with the static field applied along the c axis. The vortex dynamics are described by an 
overdamped Langevin equation with a harmonic elastic free energy. We find that the variation of the 
relaxation with temperature, average magnetic field, and local field is consistent with experiments; 
however, the calculated time dependence is different from what has been measured and the value 
of the rates calculated is roughly two orders of magnitude slower than what is observed. Combined 
with the strong experimental evidence pointing to vortex motion as the dominant mechanism for T2 
relaxation, these results call into question a prior conclusion that vortex motion is not significant in 
Ti measurements in the vortex state. 
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I. INTRODUCTION 

The motivation for this work is to test a description 
of vortex motion in the cuprates. This description, ex- 
plained below, has been used^ to dismiss the importance 
of vortex dynamics in Cu Ti relaxation in YBa2Cu307. 
Here it is used to calculate T2 relaxation, widely believed 
to be dominated by vortex vibrations. The incomplete 
agreement between our results and those of experiments 
has implications for our understanding of vortex dynam- 
ics and for conclusions which have been drawn from Ti 
data. 

Many properties of the superconducting state of the 
cuprates are remarkably conventional, while the nor- 
mal state remains very poorly understood. The vor- 
tex state in these materials, with coexisting supercon- 
ducting and normal regions, offers an attractive window 
through which to examine the normal state as well as 
excitations of the superconducting state. NMR has a 
long history of contributions to our understanding of 
superconductivityi^ A major strength of NMR is that 
many different types of measurements can be performed 
on the same sample with the same apparatus, providing 
extensive consistency checks of theory. 

NMR Ti measurements in the vortex state of high 
temperature superconductors have been used exten- 
sively to study quasiparticle excitations and vortex core 
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A key advance was the realization- that 
the unique field variation in the vortex lattice state cre- 
ated a strong correspondence between local field and 
position relative to a vortex, allowing position depen- 
dent measurements of relaxation rates. However, quasi- 
particles are not the only relaxation mechanism with 
this position dependence; there are also vortex vibra- 
tions. An early position averaged calculation^ suggested 
that vortex vibrations were the dominant contributor 
to TI Ti relaxation in Tl2Ba2Cu06. However, another 
calculation^ indicated that the more three-dimensional 
nature of YBa2Cu307 and correspondingly greater stiff- 
ness of the vortex lattice meant vortex vibrations made 
a negligible contribution to Cu Ti in this material. Both 
calculations were based on a particular model of vortex 
motion. Transverse NMR relaxation provides a forum for 



evaluating the validity of this model. 

T2 measurements in the vortex state have been 
used to study vortex dynamics near the vor- 
tex melting transition, complementing lineshape 
measurements liSiiiiiSiiiiiii^ A diffusive model of the 
dynamics has been suggested in this regime»iSi Our mo- 
tivation has been to contribute to the understanding of 
low-energy electronic excitations, which are most easily 
distinguished at very low temperatures. We therefore 
focus on temperatures well below the melting transition. 
Even in this regime, oxygen T2 relaxation in the mixed 
state of YBa2Cu307 is thought to be dominated by 
vortex vibrations for several reasons. First, the rates 
observed at the planar and apical sites are very similai^ 
as are those for yttriumifl. The dominant T2 mechanism 
for the planar oxygen in the normal state is believed 
to be dipole coupling with the Cu spinsi^. However, 
the apical oxygen and the yttrium have a very different 
crystallographic environment from the planar oxygens 
and would therefore not be expected to experience the 
same rate of dipole-induced relaxation, whereas the 
vortex-motion-induced relaxation could be the same at 
these different locations due to the continuity of the vor- 
tices. Second, the Cu dipole coupling mechanism gives 
rise to a Gaussian decay of the spin echo, whereas the 
echo decays seen in the vortex state are exponential j^iii 
Finally, the rate of relaxation varies as a function of 
local field, being faster at higher local fields and hence 
closer to the vortex coieam, strongly suggesting a vortex 
related mechanism. 

To our knowledge, no analytic calculations of T2 in 
the vortex state have been done. Ryu and StroudJi per- 
formed numerical Langevin dynamics simulations of the 
vortex state from which time scales relevant to NMR were 
extracted. Experimental results have generally been ex- 
amined assuming the field-field correlations can be char- 
acterized by a single time constant thought to be propor- 
tional to the spectral density of the field fluctuations at 
some low frequency. The accuracy of this approach has 
been shown to be poor even in fairly simple models^ 

We present here (Section^ a calculation of the rate 
of transverse relaxation due to vortex motion assuming 
overdamped harmonic vibrations. The contributions of 
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all modes are included, not just those at a single fre- 
quency. We explore (Section the time dependence 
and the local field (position) dependence of the echo 
decay as well as the temperature and field dependence 
of the characteristic decay time. The variation of our 
results with position, temperature and magnetic field 
agrees with available data. However, the time depen- 
dence predicted by our calculations is different from that 
observed, and the absolute value of the rates we calculate 
are at least two orders of magnitude slower than those 
seen in experiment. We discuss (Section llV|) possible ex- 
planations for this discrepancy and implications. A clear 
understanding of vortex dynamics is relevant not only to 
the analysis of Ti measurements as discussed above but 
also to many fundamental questions such as the nature 
of vortex core excitations and the significance of Nernst 
effect measurementaiSi as well as applied issues such as 
efficient enhancement of critical currents. 



II. CALCULATION 

In a spin echo measurement the sample sits in a 
static field Ho — HqZ and two pulses of a time varying 
field Hi are applied. Hi is in the xy plane (perpendicu- 
lar to Ho) and oscillates at the resonance frequency of a 
particular nuclear transition. The first pulse rotates the 
magnetization associated with that transition into the xy 
plane. The spins which make up the magnetization are 
initially all aligned, but due to slight variations in local 
field which cause slight variations in precession rate, the 
spins quickly fan out, causing the net magnetization to 
vanish. The second Hi pulse, applied at a time r after 
the first pulse, rotates the spins 180° about an axis j ro- 
tating at the average rate of precession. Slower spins are 
now ahead and faster ones behind, and at a time 2r after 
the initial Hi pulse the spins realign, causing a peak in 
the net magnetization. The transverse relaxation time, 
T2 , is the time scale for the decay of this spin echo with 
delay time r and represents dephasing due specifically 
to time-dependent local field variations. The effect of 
static field variations is removed in the spin echo pro- 
cess. Longitudinal (Ti) relaxation contributes in princi- 
ple to the spin echo decay; however, Ti is generally much 
greater than T2 and hence this effect can be neglected. 
(Ti - lOOOTs for planar oxygen in YBC0.3) 

We calculate here the magnitude of the net magnetiza- 
tion in the plane perpendicular to Ho at the time of the 
spin echo, 2t. Our calculation is made possible by assum- 
ing that the probability distribution of the angle through 
which a given spin has precessed, (f), is Gaussian, If this 
is the case, 



M(2r) = Afoe-<'^'(2r))/2 



(1) 



where Mo is the initial magnetization and {(fP'{2T)) is 
the thermal average of the square of the precession an- 
gle at the time 2 t. Because the instantaneous precession 



frequencies are continuous in time, this Gaussian approx- 
imation is expected to be fairly accurate. is simply 
proportional to the time integral of the local field. Tak- 
ing the effect of the second Hi pulse into account (and 
assuming its duration is negligible), the second moment 
of the precession angle may be written in terms of the 
field-field correlation function: 



2tn ^ 



dt dt'{K{t)h,{t')). (2) 



where 7^ is the nuclear gyromagnetic ratio. 

The rest of the calculation relies on a description of 
the dynamics of the local field for which we have made 
the following fairly standardiiSSiS^ assumptions. First, 
because all measurements to date have been done on 
YBCO which has moderate anisotropy, we use a three- 
dimensional anisotropic continuous model rather than 
a Lawrence-Doniach description of the superconductor. 
We keep c-axis versus plane anisotropy, but neglect ab 
anisotropy. We consider a static field [Ho) applied par- 
allel to the c axis and assume the equilibrium positions 
of the vortices form a hexagonal lattice. The vortex field 
profile is calculated from London's equation with a Gaus- 
sian short-distance cutoff. 
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ue(k, t)k± + ut(k, t)z x fc(f) ) 



r° is the equilibrium position, and ut is the horizontal 
displacement of vortex i. ui and Ut are the longitudinal 
and transverse Fourier components of the displacement. 
Xab and ^ab are the magnetic penetration depth and the 
superconducting correlation length in the ab plane, — 
^ab/r and Ac = FAafc. 

The vortex dynamics are described by a Langevin 
equation: Each vortex segment (length d) feels an elas- 
tic restoring force given by the derivative of the vortex 
lattice free energy, a viscous force proportional to its ve- 
locity with constant of proportionality 77, and a random 
thermal noise term. 



STelastic ,dUi 



elastic 



dk 
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d^, = 
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(7) 



-et(fc±,A:.)|Mt(k,t)|n(8) 



The inertial term is neglected. Details of the vor- 
tex lattice elastic properties are described in a review 
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by E. Brandt™. From these come the displacement- 
displacement correlations. 



where ji{q±,qz) 



Br] 



-7«(q)l*-t'l 



(9) 



(10) 



where B is the average internal field and $o is the flux 
quantum, and similarly for the transverse correlations. 

In calculating the field-field correlation function, we 
sum over equivalent vortex lattice positions, i, but retain 
position dependence within the unit cell, p. 

G{r,t-t') = {K{v,t)K{v,t')) (11) 

G(p,t-0 = l-Y^l-JdzG{v,t-t') 

i 

1 NT-NT^ ^,„„- f dk 



BZ 



(27r)3 



g g' 

x{h,{g' + k,t)h,{g' + k-g,t')){U) 

Here g and g' are vortex lattice reciprocal lattice vec- 
tors, and k is limited to the first Brillouin zone. {hz{g' + 
k, t)hz{g' +yi—g, t')) can be written in terms of the longi- 
tudinal and transverse displacement-displacement corre- 
lations using the harmonic form of the elastic free energy 
to write 

^g-i(g'+k)-ri(t)g-i(g-g'-k)-r,(t')^ 

^ g-jk.r°g+jk.r°g-i([(g'+k).u,(z.,t)]^) 

xe-^([(g-g'-k)-u,(^,.t')]'> 

Xe-<[(s'+k)-Ui(z,,t)][(g-g'-k).Uj(zj ,*')]) ("23) 

The Debye- Waller factors may be simplified by neglecting 
the hexagonal anisotropy of the vortex motions: 



g-i([(g'+k).u,(^.,t)]^> 



-i|g'+k|^(u^) 



(14) 



where (m^) is the mean square deviation of a vor- 
tex from its equilibrium position. We expand 

e-{[(g'+k)-Ui(zi,t)][(g-g'-k)-Uj(2j,t')]> ^ 1 _ ^[(g' _|_ . 

Ui{zi, t)][(g — g' — k) ■ Uj{zj,t')]). The zeroth-order term 
is independent of time and hence gives zero when the 
time integral is performed, leaving only the linear term. 
This appears to be analogous to the one-phonon approx- 
imation from neutron scatteringr24 However, because the 
physical meaning-and magnitude-of the quantities in- 
volved are not directly analogous, we have verified the 
approximation numericallyiS^ 

Using displacement-displacement correlations derived 
from the Langevin dynamics and performing the time 
integral, we obtain the following expression for the second 
moment of the phase distibution: 
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(-3 + 2jtT + Ae-^*'' - e-27tr-) " 



(15) 



Unlike in more strongly layered materials,^ this can- 
not be significantly simplified by physically reasonable 
approximations. We have therefore calculated the sums 
numerically. An adaptive step size approach was neces- 
sary in order to obtain reliable results in a reasonable 
amount of time given the very fast variation at small 
k values and slow variation at higher momenta with a 
time dependent crossover between the two regimes. For 
ease of calculation, a rectangular unit cell was used in k- 
space. The symmetry of our results is unaffected by this 
approximation, implying that very high k values, where 
our harmonic approximation is least justified, do not play 
a key role. Each sum on g and g' includes all recipro- 
cal lattice vectors which make a significant contribution. 
This means magnitudes of g up to roughly 27r/^afc, cor- 
responding to the short-distance cutoff. Values of g' be- 
tween g and the origin dominate, and concentric rings 
sufficient to obtain convergence are included. 

The model includes two externally controlled param- 
eters, temperature and field, for which we used values 
consistent with available experiments: temperatures be- 
tween 10 and 50 K and average internal fields between 5 
and 25 T. There are five material parameters required, 
for which we used the following values. The gyromag- 



netic ratio of ^^O, 
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3.627 X 10^ 1/sT. For the in 



plane magnetic penetration depth and superconducting 
correlation length of YBCO we used 1600 A and 16 A 
respectively, with an anisotropy ratio F of 5. The vortex 
viscosity coefficient, 77, is inferred from complex surface 
impedance measurements to be of order 10~^ kg/msi^ 
From these parameters, the mean square deviation of the 
vortices, (u^) may be calculatedi^ As a check on our 
model and parameters, we calculated (u^) at tempera- 
tures and fields where vortex melting is observed and 
found values consistent with the Lindemann criterion. 



III. RESULTS 

Two main factors affect the rate of relaxation. First, 
the larger the gradients in the local field, the faster de- 
phasing will occur, simply because it is the field varia- 
tion which gives rise to a range of precession frequencies. 
Second, the closer the match between the characteristic 
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FIG. 1: Calculated echo decay at B =10 T and T =40 K 
using model and parameters described in text. Gaussian and 
single exponential lines are shown for comparison. 



FIG. 2: Calculated transverse relaxation rate as a function of 
average mangetic field, B, at T = 20 K. 



frequencies of the vortices, 7^ and 74, and NMR time 
scales the shorter T2. If the vortices move very slowly 
the field gradient is effectively static and hence does not 
cause decay of the spin echo. If the vortices move very 
quickly, each nucleus experiences effectively an average 
field, and again there is no time dependence. Our calcu- 
lation suggests that the motions are on the fast side of 
the maximum between these two limits. Contributions 
from the transverse modes are roughly two orders of mag- 
nitude greater than those from the longitudinal modes. 
The transverse modes are softer and hence correspond to 
lower frequencies. Likewise, although all fc values con- 
tribute, it is the small k (long wavelength) modes which 
dominate at long times. 

Figure ^ shows the spin echo height as a function of 
time 2t following the initial Hi pulse. The results shown 
were calculated at 40 K and an average internal field of 
10 T and are representative of the qualitative behavior 
seen throughout the temperature and field range we ex- 
amined. The decay is neither Gaussian, as in the Cu 
T\ mechanism^^, nor simple exponential, as seen in mea- 
surements in the vortex state^^. 

We define to be the time at which the echo decays to 
1/e its initial value. Figure |21 plots I/T2 as a function of 
average internal field at a temperature of 20 K. The rate 
decreases with increasing magnetic field. This represents 
a combination of two effects: First, the field distribution 
becomes somewhat smoother at higher fields causing less 
local field variation for the same magnitude of vortex mo- 
tion. Second, the vortex lattice is stiffer at higher fields 
causing a greater mismatch between the characteristic 
vortex frequencies and the NMR time scale. A decline in 
the rate of exponential echo decay with increasing mag- 
netic field was observed by the Bachman, et alii. 

Figure 131 plots I/T2 as a function of temperature at an 
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FIG. 3: Calculated relaxation rate as a function of tempera- 
ture at an average field i? = 10 T. 



average internal field of 10 T. The rate increases essen- 
tially linearly with temperature. This represents greater 
occupation of vibrational modes at higher temperatures, 
i.e. larger amplitude motions. Curro, et al"^ show a linear 
increase in the transverse relaxation rate-interrupted at 
30 K by a peak which appears to be field independent 
and hence not vortex related. 

All the results presented above represent averages over 
the whole sample. Figure 21 shows the height of the spin 
echo as a function of position along the line between two 
vortices at 20 K, an average internal field of 10 T and 
a series of delay times. At the center of each vortex 
and at the saddle point between them, the field gradient 
goes to zero and the relaxation is vanishingly small. At 
a distance just over i^^h from the center of each vortex. 
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than circular symmetry of the vortex unit cell. The range 
of rates at a given field is, however, small. 

Curro, et al^ measured I/T2 as a function of position 
in their significantly broadened resonance line. Their re- 
sults show a rise by a factor in the range 1.5 to 2 be- 
tween the lowest and highest local fields. The broaden- 
ing of their resonance line implies an averaging of adja- 
cent rates, which is consistent with their having measured 
nonzero I/T2 at the low field end. Furthermore, the loss 
of signal at high fields would have made it difficult to 
see the downturn in the high field rates, although some 
temperatures show leveling off at high fields. 



IV. DISCUSSION 



FIG. 4: Calculated echo magnitude as a function of position 
between two neighboring vortices at B —10 T, T —20 K and 
a series of times, 2r, shown in key. 
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FIG. 5: Points show calculated transverse relaxation rate as a 
function of local field. Filled circles correspond to positions on 
a line between neighboring vortices. Open squares correspond 
to positions along a line passing through a vortex and an 
adjacent field minimum. The star is a point half way between 
the field minimum and the saddle point. The dashed line 
sketches the theoretical vortex lattice resonance lineshape. 



where the field gradient is very large, the relaxation is 
fastest. 

NMR measurements cannot look directly at this posi- 
tion dependence, but a strong correlation can be made 
between position and local field. Figure shows I/T2 
as a function of local field, again at 20 K and an aver- 
age internal field of 10 T. The dashed line sketches the 
theoretical resonance lineshape. Note that there can be 
multiple rates corresponding to a single local field. This 
is both because the rate depends on the field gradient and 
not the field itself and because of the hexagonal rather 



The trends we calculate are consistent with the tem- 
perature, average field and local field (i.e. position) de- 
pendence seen in experiments. However, the time de- 
pendence of the echo decay is different from that which 
is observed, and there is a major mismatch between the 
rates we calculate (of order 10 1/s) and the rates seen 
in experiments (of order 1000 1/s) . A mismatch of this 
magnitude was noted earlier by experimentalists^iifl com- 
paring single mode models with numerical calculationsii. 
We have shown that this discrepancy is inherent to the 
overdamped vortex vibration model and not simply cre- 
ated by the approximations involved in the single-mode 
comparison. Possible explanations for the discrepancy 
include (i) that the physical parameters we use are inac- 
curate, (ii) that our model of vortex dynamics is inappro- 
priate, and (iii) that vortices are not in fact the dominant 
mechanism for transverse relaxation. The third possibil- 
ity seems especially unlikely for all the reasons given in 
the introduction. 

Regarding the first possibility, we have studied the sen- 
sitivity of our results to the four material parameters F, 
£.ab, Aab, and 77. The anisotropy enters primarily through 
the vortex lattice elastic constants^^, and the dependence 
is extremely weak: doubling F changes the rate by less 
than 5%. The primary role of the correlation length is in 
determining the field distribution near the core. Smaller 
values of mean larger field gradients and hence faster 
relaxation. However, cutting the correlation length from 
16 A to 10 A reduces T2 by roughly 20%. The magnetic 
penetration depth influences both the field gradient and 
the stiffness of the lattice. A larger \ab produces less 
variation in the local field and a softer vortex lattice. Of 
these two competing effects, the first is stronger. Reduc- 
ing \ab from 1600 A to 1000 A reduces T2 by almost 50%. 
Finally, the vortex viscosity influences the characteristic 
frequencies. As discussed above, both very high frequen- 
cies (small 77) and very low frequencies (large 77) produce 
slow relaxation. An increase in rj from 10~^ kg/ms to 
10~^ kg/ms slows the vortices and reduces T2 by about 
50%. Even with all of these changes combined, our calcu- 
lated rate is at best still more than an order of magnitude 
slower than that observed. 
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Regarding the possibility that the vortex dynamics are 
not captured by our model, the overdamped Langevin dy- 
namics picture is widely usedpili^SiSi Undamped motion 
was explored in the context of organic superconductors,^^ 
but would produce even higher characteristic frequencies. 
Another possibility might be a slower motion superim- 
posed on the faster motion we've described, for exam- 
ple flux creep of vortex bundles. However, stimulated 
echo measurementsifi suggest that the vortex motions are 
small compared to their lattice spacing. As for the influ- 
ence of disorder, the average intervortex spacing is not 
dramatically changed^^ by disorder and our results ap- 
pear not to be very sensitive to short wavelength distur- 
bances. 

A key feature of our implementation of this model, 
however, appears to be the assumption of a continuum 
of vibrational modes, corresponding to an infinite sys- 
tem size. The influence this has on the time dependence 
can be seen as follows. For each mode, when 7r << 1 
{(fP') oc and when 7T >> 1 {(jp) (x r. Because in our 
model there is a continuum of modes down to zero en- 
ergy, one never reaches the long time limit of all modes 
and therefore there is always curvature. Available vibra- 
tional modes will be restricted in a finite size system, 
such as the powder samples used in NMR, and also in a 
system in which vortex pinning is significant. A very sim- 
ple model in which short wavelength modes are simply 
removed does indeed produce an exponential echo decay 



as seen in experiments. However, a more careful explo- 
ration of the effect of a discrete spectrum on our results 
is still in progress. 

In conclusion, we have calculated the time, tempera- 
ture, average field, and local field (position) dependence 
of the rate of transverse NMR relaxation arising from 
vortex motion using an overdamped Langevin dynamics 
model and harmonic elastic free energy. The functional 
dependence of our results on temperature, field and posi- 
tion are consistent with available experiments, while the 
time dependence and the rate itself are not. A key im- 
plication of this result is that vortex motion cannot nec- 
essarily be neglected as a mechanism for Ti relaxation 
in the vortex state, influencing conclusions which have 
been drawn from these measurements on the nature of 
the electronic states in the presence of vortices. 
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